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^ ! Abstract 

Production of /q(1710), a theoretical endeavor of pure scalar glueball state, is studied in de- 

tail from exclusive rare B decay within the framework of perturbative QCD. The branching 

fraction for B ± -> K^f (1710) -► K*±(KK) is estimated to be about 8 x 10~ 6 , while for 

\ B^ — > K ^ fo (1710) — > K^(KK) it is smaller by roughly an order of magnitude. With the ac- 

(■ — | cumulation of almost 1 billion BB pairs from the BaBar and Belle experiments to date, hunting 

\Q ' 

■ for a scalar glueball via these rare decay modes should be attainable. 
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From the modern point of view, properties of pseudoscalar mesons can be understood as 
Nambu-Goldstone bosons due to the spontaneous symmetry breaking of chiral symmetry. 
Their low energy dynamics can be described by the chiral lagrangian. On the other hand, 
scalar mesons are not governed by any low energy symmetry like chiral symmetry and thus 
they can not take advantage of the power of a symmetry. Indeed, their SU(3) classification, 



the quark content o 



or sea 



ar mesons [l]. Moreover, possible mixings of the qq states with a pure glueball state 



their composition, as well as their spectroscopy are not well understood 



must be taken into consideration. 
Recent quenched lattice simulation |9|] predicted the lowest pure glueball state has a mass 
equals 1710 ± 50 ± 80 MeV and J PC = ++ . The first error is statistical while the second is 
due to approximate anisotropy of the lattice. This suggests that before mixing, a glueball 
mass should be closed to 1710 MeV, instead of the earlier lattice result of 1500 MeV 
This makes / (1710) a strong candidate for a lowest pure glueball state as advocated in 
[h]] based on argument of chiral suppression in / (1710) decays into pair of pseudoscalar 
mesons. The next two pure glueball states predicted by the quenched approximation 9] have 
masses 2390 ±30 ±120 MeV and 2560 ±35 ±120 MeV with J PC = 2 ++ and 0"+ respectively. 
Mixings between the nearby three isosinglet scalars / (1370), /o(1500), and / (1710) and the 
isovector scalar ao(1450) have been studied in detail in [2| with the following main result: In 
the SU(3) symmetry limit, / (1500) is a pure SU(3) octet and degenerate with the isovector 
scalar meson ao(1450), whereas /q(1370) is mainly a SU(3) singlet with a small mixing with 
/o(1710) which is composed predominantly by a scalar glueball. 
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Important production mechanism of glueballs is the decay of heavy quarkonium 
In fact, the observed enhancement of the mode J /if) — > / (1710)o; relative to / (1710)</> 
and the copious production of /o(1710) in the radiative J/tf> decays are strong indication 



that / (1710) is mainly composed of gluebal 



a. 



Another interesting mechanism is the 



direct production from e + e — > 7* — > GjH jl4|. where Gj stands for a glueball state of spin 
J = or 2 and H denotes a J /if) or T. Recently, glueball production from inclusive rare B 



decay [15] has also been studied. Ironically, scalar glueball state has never been observed in 
the gluon-rich channels of J/if){lS) decays or 77 collisions 1 . 

In this article, we will study glueball production via exclusive B decay using perturbative 



1 For a summary of the non-qq candidates from the Particle Data Group, see p949 of Rcf. 16]. 
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QCD (PQCD). Firstly, we will ignore mixing effects and treat /o(1710) as a pure scalar 
glueball suggested by the quenched lattice data. At the end of the paper, we will demonstrate 
the mixing effects are minuscule. At quark level, the effective Hamiltonian for the decay 



b — ► sqq can be written as 

G. 



H, 



eil' 



— V V 



q=u,c 
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8=3 



(1) 



where V q = V*V q b denotes the product of Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements and the operators Oi-Oio are defined as 



= {s a q/3)v-A{qpb a )v-A , 



Of 



o 3 
o 5 



[s a b. 



a V-A 



IV-A , 



[s n b. 



a u a)V-A 



q 



K s a q a )v-A{(lpbj})v-A , 
4 = {s a bp)v-A^2,{qpq a )v-A , 



[s n b 



"" 11 V-A /"] 

q 



a)V+A , 



7 = -{s a b a )v- A ^2 e q(Wn)v+A , 8 = -{s a b p ) v „ A ^ e q {q p q a \ 

q q 
3 3 

09 = -(s Q &a)i/- J 4^e g (g /3 g / g)y-A , Oio = 7T {sah) v-A 7^ e g {qpqo 



V+A , 



IV-A 



(2) 



with a and (3 being the color indices and 0i-0io the corresponding Wilson coefficients. 
In addition, the gluonic penguin vertex for b(p) — > s(p')g*(q) with next-to- leading QCD 



corrections is given by 



Gr (Js V t *Vtbs(p') [AF L (g 2 7 Al - q 1 * i)L - zm b F 2 cr^q u R]T a b(p), 



(3) 



where g s is the strong coupling constant, mj is the 6-quark mass, T a is the generator for 
the color group, and L(R) = (1 =F7s)/2, AFj = An(Cf(q, fj.) + Cf(q,fi))/aJu) and F 2 = 
-2GgV). Explicit formulas for 0f , 0f , and 0g can be found in Ref. Q. Since the 
ground state scalar glueball is composed of two gluons, the effective interaction between a 
scalar glueball and gluons can be written as [10 ] 



(4) 



where /o stands for an unknown effective coupling constant, Go denotes the scalar glueball 
field, and is the gluon field strength tensor. With these 4-quarks operators 0\ — O w 
(12]) and the two effective couplings (j3j) and (jl]), we can embark upon the computation of the 



decay rates for B — ► K^*'Gq using PQCD. The flavor diagrams for B — > K^*>Gq decays are 
displayed in Fig. [U Fig. (la) denotes contribution from the 4-quarks operators Oi-io given 
in Eq.(121), whereas Fig. (lb) involves the gluonic penguin vertex contribution of Eq.Q. Both 
diagrams are of the same order in a s . In the heavy quark limit, the production of light meson 
is supposed to respect color transparency Q, i.e, final state interactions are subleadine 
effects and negligible. We will work under this assumption in what follows. Moreover, 
diagrams like FigfS] that are of higher order in a s will be ignored. 





FIG. 1: Flavor diagrams for the B -> K^G . 





FIG. 2: Other flavor diagrams for the B — ► K^*'Gq at higher order in a s 



To deal with the transition matrix elements for exclusive B decays, we employ PQCD 



21 



221 ] factorization formalism to estimate the hadronic effects. By the factorization the- 
orem, the transition amplitude can be written as the convolution of hadronic distribution 
amplitudes and the hard amplitude of the valence quarks, in which the distribution am- 
plitudes absorb the infrared divergences and represent the effects of nonperturbative QCD. 
As usual, the hard amplitudes can be calculated perturbatively by following the Feynman 



rules. The nonperturbative objects can be described by the nonlocal matrix elements and 
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are expressed as [23l. \24. 125] 

(VZ e-* kz (Q\bp(Q)q a (z)\B(p B )) = --L=[(fi B + m B ) lb Ucj> B {k) 



(2tt) 4 x 1 pv " v^at: 

il ' Z e-^ z (K(p K )\qp(z)s a (0)\0) = -^={[75 ArUM*) + HUm° K f K {x) 



+m° K [^ 5 (A+ A- ~ 1)Wk(z)} 
f/ '* e -^— (ir*(p^,e i )|^(z)s ce (0)|0) = -i={ m ^[^] Q ^( x ) 



(2tt) 4 n ^ 

+ [A i>K*]ap4^ K *{ x ) + rriK^lla^^x)} ,(5) 



for 5, K, and mesons, respectively, where N c is the number of color, n± are two light-like 
vectors satisfying n+ -n_ = 2, and is the longitudinal polarization vector of K* . <j) B (x, b), 
the distribution amplitude of B meson, is constructed as follows [251 ] 

(j) B (x,b) = J dk + d 2 k ± ^^(f)B{k) (6) 

with x = k~ /pg. (ft K (*)(x) and 4> p ^ (t) {x) are the twist- 2 and 3 distribution amplitudes of 
mesons with the argument x stands for the momentum fraction. Finally, m B and m K (*) 
are the masses for the B and K^*> with m° K = m 2 K /(m s + m q ) where m q and m s denote 
the light quark masses. The meson distribution amplitudes are subjected to the following 
normalization conditions 

fB(K(*)) f 1 A ^ AP fxl) 



dx<f) B(KM) {x) = , J dx<j) P KM (x) = ^-J=, ^ dx<f> a KW {x) = (7) 

where <p B (x) = 4> B (x, 0) and fsucW) an d fjc(*) are ^ ne decay constants. We do not introduce 
transverse momenta for the light mesons K and K* here which we will justify later when 
we discuss the end-point singularities of the decay amplitudes. 

In the light-cone coordinate system, we can pick the two light-like vectors to be n + = 
(1, 0, 0j_) and n_ = (0, 1, 0j_), and the momenta of the B and K mesons can be written as 

p B = ^|(1,1,0jl), px = ^(1-^X1,0,01), (8) 

with tq = mG /m B . For the vector meson K*, we take 



PK* = ^(l-r 2 Go ,r 2 K *,0 ± ), e L = -J—(l-r 2 Go ,-r 2 K „0 ± ), (9) 
V2 V%r K * 



with rx* = mK*/fn B in which the physical condition el-Pk* = is satisfied for massive vector 
particle. The momenta of the spectator quarks with their transverse momenta included are 



given by 

h = (o, hxj , k 2 = y^i 1 - r G ) x 2,0, k 2 ^j . (10) 

With these light-cone coordinates and distribution amplitudes defined, we can study the 
transition matrix elements for B — > MGq (M = K, K*). We first analyze Fig. [T](a). Within 
the PQCD approach, we find that Fig. Q^a) is directly proportional to x\. Since x\ is the 
momentum fraction of the valence quark inside the B meson and its value is expected to be 
x\ ~ kjvriB <C 1 with A = ms — nib- Comparing to x 2 ~ 0(1) (Fig. Hfb)), its contribution 
belongs to higher power in heavy quark expansion. As an illustration, we can use the 
operator O4 in Eq. (j2J) to demonstrate this effect. Thus, one finds 

M 04 oc ^ o9 ^p F y/N~ c f K m% [ ^1^2777^777^5(1 y-\ (2 - x 2 )x 2 xi0 B (xi, A?i_l) 

v2 J {^) {^) V m B J 
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(m 2 Go - m\x x ) - \k xl \ 2 m 2 B x 1 x 2 - \k 1± - k 2 ±\ 2 m 2 Go {l - x 2 ) - \k 21 _\ 2 I^LI 2 ' 

(11) 

where Cf = {N 2 — 1)/(2N C ). It has been shown in [25( that under Sudakov suppression 
arising from k± and threshold resummations, the average transverse momenta of valence 
quarks are (k±) ~ 1.5 GeV and the end point singularities at x\ y2 — > in Eq. tflTT) can be 
effectively removed. With an explicit factor of x± appearing in the numerator of Eq. ffTTl) . this 
contribution is regarded as a higher power effect in l/ms and therefore can be neglected. We 
note that this situation is quite similar to the flavor singlet mechanism to the B — > 77' form 



factor 26|]. According to the PQCD analysis in Ref. 27|, contribution from the possible 
gluonic component inside 77' to the B — > 77' form factor also has similar behavior. Its 
numerical value is two orders of magnitude smaller than the B — > tt form factor. Similarly, 
other operators Oi_3 and Os_io give the same behavior. Therefore, to the leading power in 
Aqcd/^b, the contributions from Fig. [T](a) can be neglected. We will concentrate on the 
contribution of Fig. QJb) in what follows. 

By using the introduced nonlocal matrix elements for mesons and the light-cone coor- 
dinates given above, the transition matrix element for B — ► MG (M = K, K*) can be 
obtained from Fig. HJb) as 

A M = ^kv ts V t * b M M (12) 



with the decay amplitude function M.m given by 



M.M = — —JqCf / bdb / dx\ I dx 2 <pB(^i,b) 
K Jo Jo Jo 

xx 2 (e^^> M (x 2 ) + ef}4> p M (x 2 ) + (j) a M {x 2 )\ E(t)h(x u x 2 , b) 

e® = m(t)(l - r 2 Go )[l + 2r 2 Go + 2(1 - r Go )x 2 ] - 3r 6 (l - r 2 Go )F 2 (t) , 
eg } = 3r K [-2AF 1 (t)(r 2 Go + (1 - r Go )x 2 ) + r b F 2 (t)(l + r 2 Go + (1 - t 2 Gq )x 2 )] 
e% ] = r b r K (l - r 2 Go )(l - x 2 )F 2 (t) , 

for the pseudoscalar K, and 



(13) 



'if* 



(1) (2) _ ^fl n _ 2 N (3) (3) 



T-Jf* (2) 
e _fsT > 



(14) 



(15) 



for the vector meson X*. Here we have introduced the dimensionless variables r& = m^/mB, 
rx = m° K /mB, and r^* = m^/ms- The hard function h{x\,x 2 ,b) in Eq. (|T3|) is given by 

1 



h(x 1 ,x 2 ,b) 



X 12 +Y 1 



12 



K 



m 2 B Y 12 b 



m 2 B X 12 b 



(16) 



with X\ 2 = (1 — Xx)[r Go + (1 — r Go )x 2 ] and Y12 = (1 — r Go )xix 2 . The evolution factor E(t) 
in Eq. ljTBl is defined by 



E{t) = a s (t)e 



-S B (t)-S K (t) 



(17) 



where exp(— Sb(k)) is the Sudakov exponents that resummed large logarithmic corrections 
to the B(K) meson wave functions 28|, [29|. Their explicit forms are given by 



5 /"* <i/i 

S B (t) = s{x x p%,b) + - / — 7(a a (/2)) 



S K {t) = s(x 2P + , b) + s((l - x 2 )p + , 6) + 2 



— 7(a,0x)) 
1/6 Z 1 



(18) 



where 7(a s (//)) is the anomalous dimension. To leading order in a s , j(a s (^)) equals —a s /ir. 
The function s(Q,b) in Eq. ffTBl is given by [30, [31] 

•Q 



A(a s (/x)) + B(a s (iJ,)) 



(19) 



where 



with / = 4 being the active flavor number and 7^ is the Euler constant. As mentioned 
before, x\ ~ A/m B <C 1, we have dropped all terms related to x\ in the above expressions 
for {e^}. Since r K {*) < 1, we have retained only those terms in the above formulas for {e^} 
that are at most linear in r K (*). The scale t where the strong coupling a s (t) in (|17[) . the 
Sudakov exponents in (IT5T) . and the AFi(t) and F 2 (t) in (TT4"j) are evaluated will be discussed 
later. For comparison, we also present the formula of the decay amplitude function M.m 
with k± = in Appendix lAl 

For estimating our numerical results, we take the values of theoretical parameters to 
be: f B = 190 MeV, m b = 4.4 GeV, (m B , m K , m K *, m Go ) = (5.28, 0.493, 0.892, 1.71) GeV, 
V ts V t * b = —0.041. For the B meson distribution amplitude, we use [28]] 

Mxub)=N B xi(l- Xi yexp Y-;U& 2 (21) 

with N B = 111.2 GeV and uj b = 0.38 GeV. For the distribution amplitudes of the light 
pseudoscalar K and vector mesons K*, we refer to their results derived by the light-cone 



QCD sum rules in [32|, |33j, [34J . Their explicit expressions and relevant values of parameters 
are collected in the Appendix [B] for convenience. 

According to the results of light-cone QCD sum rules, at small the behavior of twist-2 
distribution amplitude obeys the asymptotic form <p K ^(x 2 ) oc x 2 (l — x 2 ), whilst those of 
twist-3 distribution amplitudes approach a constant 0^ t) (x 2 ) oc const. Consequently, at 
small x 2 , the decay amplitude function contributed by the twist-2 distribution amplitudes 
of behaves like 



tw2 x 2 (j) B (xi)4) K (*)(x 2 ) x 2 x\{\ - xi) 2 x 2 (l - x 2 ) _ xi(l - xi) 2 a;2(l - x 2 ) 
KM k 2 q 2 x lX2 (rl o + (1 - r 2 Go )x 2 ) (r Go + (1 - t 2 Gq )x 2 ) 



Obviously, even if one sets r Go to be zero, the effects from twist-2 distribution amplitudes 
of are well-defined at the end point x 2 — > 0. Similarly, the contribution from twist-3 
distribution amplitudes to the decay amplitude function at small x 2 behaves like 

tw3 x 2 B (xi)0^ ) (x 2 ) x 2 x 2 (l- Xl ) 2 _ x^l-x,) 2 

kW x lX2 (r 2 GQ + (l-r 2 Go )x 2 ) (r 2 Go + (1 - r|> 2 ) " 1 } 

Whence r G() 0, one will suffer logarithmic divergences from the twist-3 distribution ampli- 
tudes. In practice, r Go ~ 0.32, the divergence will not occur. This implies that the influence 
of k± can only be mild. As a common practice, we do not introduce transverse momenta for 
the valence quarks to suppress large effects from end point singularities. 
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Since the Wilson coefficients are \i scale dependence, for smearing its dependence, we 



include the values of Wilson coefficients with the next-to-leading QCD corrections [19]. 
However, even so, the Gf^ %g are still slightly /i-dependence. Due to this reason, determina- 
tion of the scale of exchanged hard gluons in Fig. [TJ is also one of the origins of theoretical 
uncertainties. For the gluon that attached to the penguin vertex b — > sg*, it carries a typical 
momentum of \fq^ = nis ((1 — Xi)(rQ + (1 — Tgo) X2 ) • When x% ~ A/uif, and X2 is 0(1), 
say x 2 = 0.5, we get y/q* ~ 3.9 GeV. However, the gluon attached to the spectator quark 
carries roughly a typical momentum of y/—k 2 = ttib ((1 — r^Jx^) ~ 1.4 GeV. We note 
that a suitable range of £2 in PQCD is often taken as ~ 0.3 — 0.7. For defmiteness, we take 



the democratic average value t = (\/(f + \ / —k 2 )/2 as the hard scale, in which the allowed 
value is within the range t pa 2.45 ± 0.45 GeV. This justifies somewhat the validity of the 
PQCD approach and we will take this range of t as our theoretical uncertainties. For illus- 
tration, we present the involving Wilson coefficients at different values of fi scale in Table [H 



TABLE I: The involving Wilson coefficients at various values of /i scale. 



Wilson coefficient 


fi = 2.1 GeV 


2.5 GeV 


3.0 GeV 


cf 


-(6.17 + 0.78?) x 1(T 2 


-(5.80 + 0.89f) x 10~ 2 


-(5.48 + 0.89i) x lO" 2 


/-feff 
°6 


-(7.69 + 0.78i) x 1(T 2 


-(7.19 + 0.89i) x 10- 2 


-(6.77 + 0.89?) x 10~ 2 


/-feff 
°8 9 


-0.170 


-0.165 


-0.161 



Effective interactions between a scalar glueball G and the pseudoscalars have been stud- 



ied using chiral perturbation theory 



15 



35j. By using the current experimental data [if] ] 
r to tai(/o(1710)) = 137 ± 8 MeV and BR(/ (1710) -> KK) = 0.38+gJi, this allows us to 



get an estimate of the unknown coupling /o — 0-07]to!oi8 

GeV" 1 [15]. This result of / 
should be taken as a crude estimation. For one thing, the data of the branching ratio 
BR(/ (1710) — > KK) was not used for averages, fits, limits, etc. by the PDG [if]]. Instead 
the following two ratios were used in the PDG analysis: 

r(/o(i7io) -> w) = 

nl r(/ (1710) -> KK) ' v ; 

7 r(/ (i7io) -> Kicj v ; 

Within the approach of chiral perturbation theory [l5|, it would be difficult to accommodate 
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these two ratios of Eqs. (l2l]) and (I25I) . since the leading term in the chiral Lagrangian is flavor 
blind. Here we will present another approach to estimate Jq. At quark level, the amplitude 
for G — > qq is proportional to the quark mass m q and therefore chirally suppressed. Its 



explicit form is given by 
A(G 



io| 



167TA/2m g 

qq) = -joa« — = -in 



l + {3 



U q Vq , 



(26) 



3 /3 \l-f3, 

where (3 denotes the velocity of the quark and u q (v q ) is the quark (anti-quark) spinor. It has 
been argued in [10!] that the chiral suppression of the amplitude A(Gq — > qq) oc m q persist in 
all order of a s . One may treat the coefficient of this decay amplitude as the short- distance 
coefficient of the strong decay Go ~~ ► PP where P stands for a pseudoscalar meson like ir, 
K, or r] etc, as illustrated in Fig. [31 Thus, 



(PP\H eS \G ) = -f m q YF 



pp. 



m 



Go) ' 



with, to leading order in a s 



Y 



In 



1 + P 



, .levrv^l 
Go ' 3 p"'\l-f3 



(27) 



(28) 



and F pp (mQ ) is the time-like form factor (PP \qq\0) evaluated at Q 2 = rn>Q o . For the case 



Go 




FIG. 3: Flavor diagram for the Go — * PP with P being a pseudoscalar. 

of P = rj, we include the quark-flavor mixing effect according to rj = cos(pr] q — sin0?7 s and 
rj' = sm<f)r] q + cos(pr] s with rj q = (uu + dd)/^/2, r] s = ss 3y, l37J, and <p — 41.4° [38j. Using 
Eq. (I27p . the following ratios of the partial decay rates can be obtained 

M 2 



R 



■k/K 



r(G c 



7T7T 



R 



V/K 



r(G - 
r(G 



KK) 
* r]7]) 



3 (m u + m d f \F™{mlX (l - 4m|/m|J; 
8 ml \F KK {mlX{l-Aml/ml o ) 1/2 
(1 - 4m 2 v /m 2 Go ) 1/2 



x 



1/2 

w u + m d ) cos 2 (f)F VqVq (m G() )/2 + m s sin 2 <pF VaTta (m 2 Ga )\ 2 



r(G — > KK) (i-im 2 K /m% ) 



2ml\F^{m 2 GQ ) 



(29) 
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By taking the flavor 377(3) approximation, one finds that F n7T /F KK « f 2 /f K , F^"/F KK 
P q /P K , and FW'/F™ « /*//£. With m u = m d = 10 , m s = 120, / ff = 130, f K = 160 (l6|, 
f q = 140, /, = 180 37] (all in unit of MeV), one deduces 

Rn/K = 0.006 , R v/K = 0.37 . (30) 

Identifying Go to be /o(1710), these ratios are consistent with the current experimental data 
quoted in Eqs. ([2"l"l) and (125]) . Using Eq. (1261) . the following expression of / can be obtained 

fl = ,jTf-y» f 1 - #) BR(Go - KK) , (31) 



\F KK (m 2 Go )m s Y\ 2 \ m 2 Go 

where Tg = 137 ± 8 MeV is identified as the width of /o(1710). The time-like form factor 
F KK (m G ) has been extracted in Ref. 



39( | by performing the data fitting to non-resonant 



B — > KKK decays with the following form 

F KK (Q) = V - (SF^I + 2F^ + a NR exp(-a NR Q 2 ) , 

^-{i+w)^)' 1 - (32) 

where v = {m\ - ml)/(m s - m d ), A = 0.3 GeV, = -3.26 GeV 2 , x£ ] = 5.02 GeV 4 , 
xf } = 0.47 GeV 2 , xf = 0, omr = 4.4e^/ 4 GeV, and a NR = 0.13 GeV" 2 . By using 



BR(G -> KK) = 0.38±g;?i |1Q|, the value for f is estimated to be f = 0.086±g;g^g, 
which is only slightly larger than the value obtained from the chiral Lagrangian approach. 
In passing, we note that, using light-cone distribution amplitudes, it has been argued in 



Ref. [35| that G — > irir,KK might be dominated by the 4-quark process of G — > qqqq 
which is not chirally suppressed. Using this 4-quark mechanism and PQCD factorization 
scheme, one would predict a large ratio of Rn/K ~ (/tt//x) 4 ~ 0.48. For further discussion 



of this issue, we refer our reader to Refs. 



35 



40 



4l| 



Using the matrix element defined by Eq. (|T3l with the above chosen values of parameters, 
the values of M. K {*) are given in Table HT1 for /o = 0.086 GeV -1 and three different values of 
fi scale. For comparisons, we also present the results with k± = in Table HT1 

The branching fractions for B + — > (K + , K* + )Gq decays are tabulated in Table IIITl From 
Table IIIIl we find that the branching fraction for B + — > K* + Gq is about one order of 
magnitude larger than that for B + — > K + Gq. The difference arises not only from the values 
of the decay constants fx and fx* entered in the distribution amplitudes, but also from the 
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TABLE II: Decay amplitude M M (in units of 1CT 4 ) for B+ (K + , K*+)G with and without 
k± at /o = 0.086 GeV -1 and three different choices of fi = 2.1, 2.5, and 3.0 GeV. Numbers given 
in brackets are without 



Mode 


H = 2.1 GeV 


2.5 GeV 


3.0 GeV 


K+G Q 


-3.54 - 0.42i 


-3.34 - O.Ui 


-3.22 - 0.48f 




(-3.51 -0.38i) 


(-3.28-0.41i) 


(-3.08 - 0.43i) 


K*+G 


-11.13 - 1.5H 


-12.56 - 1.51* 


-12.40 - 1.70i 




(-10.90 - 1.17i) 


(-10.18 - 1.25?) 


(-9.60 - 1.33i) 



effects of &k yk{ x ^) an d e K yk^z) m the K + Go mode, which are switched to e K 4>k*( x 2) 
and (j? K *(x2) respectively in the K* + Gq mode. We also find that the k± influence on 
B + — > K* + Gq decay is stronger than B + — > K + Gq. In addition, when [i is smaller, k± has 
lesser effects on the decay B + — > K* + Gq- 



TABLE III: Branching fractions (in units of 10~ 6 ) for B + -> (K + , K* + )G with and without k± 
at /o = 0.086 GeV -1 and three different choices of /i = 2.1, 2.5, and 3.0 GeV. Numbers given in 
brackets are without k±. 



Mode 


H = 2.1 GeV 


2.5 GeV 


3.0 GeV 


K+G 


3.05 


2.72 


2.53 




(2.99) 


(2.62) 


(2.34) 


K* + G 


29.07 


35.94 


36.06 




(26.50) 


(23.21) 


(20.69) 



The branching fractions for the decay chains B + — > K + Gq — > K + (KK)g and B + — > 
K* + Gq — > K + (KK)g are tabulated in Table IIV[ where the errors are coming from the 
experimental data of BR(/o(1710) — > KK). From Table HV] we learn that one has a better 
chance to look for the ground state of glueball through the three-body decays B — > K*KK, 
since its branching fraction could be more than a factor of 10 larger than B — ► KKK. 
Recently, BaBar had reported the following branching ratio for B ± — > (K + K~)K ± where 



the (K + K ) pair coming from the /q(1710) 42] 



BR(B ± -> (iC+K-)^^)^) = (1.7 ± 1.0 ± 0.3) x 10" 6 . (33) 
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TABLE IV: Branching fractions (in units of 10" 6 ) for B + -> (K+, K* + )(KK) Go at fjt = 2.1, 2.5, 
and 3.0 GeV. Numbers given in brackets are without k±. 



Mode 


H = 2.1 GeV 


2.5 GeV 


3.0 GeV 


K+(KK) Go 


(1-1318$ 


i.03is:« 
(i-ool8;|) 


0.9618:11 

(0-89l8f 7 ) 


K*+(KK) Go 


ii-05l|;l 

(I0.07l|;^) 


±O.OO_ 10 34 


1 3 70+ 7 - 42 
(7.861H) 



From the first and second rows in Table HVl identifying Go to be /o(1710), one can see that 
our predictions are consistent with the experimental data. 

Before we close, we want to address the issue of mixing effects. Although we have treated 
/o(1710) as a pure gluonic state, it should be interesting to consider its mixing effects with 
other qq states. To deal with the mixture of a pure glueball with the qq quarkonia state, we 
follow Ref. [2I to express the /o(1710) state as the following combination 



where \G) is the pure g. 
of the mixing schemes 



|/o(1710)) = C N \N) + C S \S) + C G \G) (34) 

ueball state, \N) = (uu + dd)/ v2, and \S) = ss. Accordingly to one 
2|, the coefficients took the following values: = 0.32, Cs = 0.18, 



and C G = 0.93. ' 
others in Refs. 0, 



"he quoted results of these coefficients are similar to those obtained by 
43I . The corresponding flavor diagrams for the decays B — > K^*\N, S) 



are shown in Fig. HI Since the distribution amplitude and decay constant for / (1710) are 
uncertain, for simplicity, we use factorization assumption to estimate the hadronic effects 
for these two-body B decays. In terms of the operators in Eq. (j2j), one can easily show that 
the contributions from diagram Fig. IH(a) and (d) are associated with the matrix element 
(N(S) I ?7 At ?|0). Since the scalar N or S is C-even while <n M g is C-odd, the contributions 
from Fig. H^a) and (d) must vanish because charge conjugation is a good quantum number 
in strong interaction. On the other hand, the contributions from Fig. IH(b) and (c) are 
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1 



1 



\J 1 




(a) 



(6) 




b 

T 

i 




(c) (d) 

FIG. 4: Flavor diagrams for the B — » K^*\N, S). (a) and (b) are from QCD and electroweak 
penguin diagrams, while (c) and (d) denote the tree contributions. 



non- vanishing and they can be derived as 
Akn 



(*F JrCn f 2 2 \ Ut t/* / « «\ t/ TA* 1 T?BN I 2 \ 



A K *s 



y/2 V% 

G p 

~V2 



2m s fsC t 



rrib + m s 



F o K ( m s) , 



^ (m B - m N ) [V ts V tb a 4 - KsKb ai] 

v/2 



m. 



(35) 
(36) 
(37) 
(38) 



for 5+ -> X+(iV, 5) and 5+ 
defined by 



m 6 - m s 

K* + (N, S) decays, respectively, where px, a±, and 



PK 

a i 



2m 2 K 



m s + m u )(m h + m u ) 



c 2 + 



C3 3 / „ Co 

C* + T7- + ^ A Cio + ^ , ■ 



C 7 



(39) 



e q is the electric charge of quark q and F with M = N, S and correspond to the 
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— I 



B — > (M, K*) form factors parametrized by 
(N( P )\h^sq\B( PB )) 
(K(p)\h^\B(p B )) 
(K*(p, e K *)\h^s\B(p B )) 



44 



45J 



m 2 B -m 2 N 



m\ - m 2 N 



n F l 



BN 1 



Pn 



m, 



m 



K 



m 



B 



m 



K 



2m K *A 



BK- 



(9 )— IT - ^ 



+ (m B + m K *)A* K * (q 2 ) (e^ - 



-A 



BE"* 



(9 2 



tub + ttik* 



Pn 



q 

P ■ q 

q z 



(40) 



Af K *(q 2 ) and A BK * (q 2 ) are two other form factors that are not relevant in our analysis. 



With 



P 



2.0 GeV, V us 
^bn 



0.22, V„ h = 3.6 x 10- 3 e-^ 3 , 



GeV [2], F BN (m 2 K ) = 0.26, F 



BN, 



m 



K* 



0.28, F BK \ 



72°, m N = 1.47 GeV, m s = 1.50 
= 0.38, A BK *(m 2 ) = 0.42 



and fs = —280 MeV 46j], one has the following estimation 

Gftti b 



Arn + Aks 
Ak*n + Ar^s 



\),\'n, (-8.50 + 1.37 i) x 10" 5 
VtsV+t (1.17 + 0.190 x 10" 4 . 



(41) 
(42) 



Comparing these values to those coming from the contribution of purely gluonic state given 
in Table [TTJ, one can conclude that the qq quarkonia contributions can be safely ignored. 

In summary, we have studied the scalar glueball production in exclusive B decays by 
using PQCD factorization approach. The typical momenta carried by the exchanged gluons 
in the process is about half of the B meson mass. One thus expects our perturbative 
results are trustworthy. According to our analysis, we find that the branching fraction for 
B + — > K + G is a few xl0~ 6 ; however, for B + — > K* + G it can be as large as 3 — 4 x 10~ 5 . 
As a result, the branching fraction for the decaying chain B + — > K^ + Gq — > K ( -*^ + (KK)g 
is ~ 0.66(7.79) x 10~ 6 . With the experimental inputs of Eqs. (l24l and (1251) . we also expect 
the branching ratios for B + — > K^ + (r]r])G and B + — > (nn)G are about 50% and 

less than 10% of B + — > (KK)g respectively. In this work, we have focused on the 

charged B mesons. Similar conclusions can be drawn for the neutral B mesons where the 
only difference is their lifetimes. Their mass difference (m^o — m B +) is merely 0.33 ± 0.28 
MeV 16j and will not affect our numerical results significantly. Thus dividing the branching 
fractions given in Table HTT1 and Table HVl by the ratio t b +/t b o = 1.071 ± 0.009 16] from 
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direct measurements, one would obtain the corresponding branching fractions for the neutral 
B meson modes. Experimentally, the mode B ± — ► (KK) f (mo)K ± has been detected at 
BaBar with a branching ratio consistent with our PQCD prediction. This work suggests 
that detection of the resonant three-body mode B — ► {KK)f (yf\_Q\K* with a predicted 
larger branching ratio can give further support of /o(1710) is a pure scalar glueball. 



APPENDIX A: DECAY AMPLITUDE FUNCTION M M WITH k ± = 

Since the mass of glueball is much larger than those of ordinary pseudoscalars, we find 
that the influence of transverse momentum on the two-body decay B — > K^Gq is not as 
large as in the case of B decays into two lighter mesons. Setting k±± and k 2 ± in the momenta 
of the spectator quarks in Eq.flSJ) to be zero, the decay amplitude function M.m given in 
Eq. flTBl) can be simplified to be 

Mm — ^ F [ dxi [ dx 2 4>b(xi) 
Km B J J 

x {^(f)M{x 2 ) + e M ) 4> p M {x 2 ) + e^M^)} a s {t)h{x 1: x 2 ) , (Al) 
with the hard function h{x\,X2) given by 

^-^- x 1 (l-, 1 )(4„ 1 + (l-r|,)x 2) - (A2) 
APPENDIX B: DISTRIBUTION AMPLITUDES FOR tfW 

In this appendix, we compile the light-cone distribution amplitudes that entered in our 
calculations. The distribution amplitudes for K, defined in Eq. (JS}, are expressed as follows 
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22| : 



<Pk{x) 
(f) p K (x) 



A 



2^/2N c 



6x(i - x) i + of cf /2 (o + c^qno 



K n 3/2, 



2^j2N c 



1 + 3p*(l + 6a?) - 9p5af + ^ /2 (£) ( ?p?af - p K _ ( + 27af 



+C 2 1/2 (0 (30 mK + lhpl4 - 3p x af ) + C 3 1/2 (0 (l^ K \ ZK - ^-p K 4 



3r]3KU} 3 KCl /2 (0 + U -{ P « + P K ) (1 - 3af + 6af ) ln(l 



+ |(p?"P A ')(l + 3af + 6af)lnx 



^f= {e [fci + fccfto + hCl /2 {i) + 6 4 C 3 3/2 (0 

— 306 3 x(l — x) + 65 ln(l — x) + 6 6 In 2 ] 

+x(l - x) [-66 2 + 56 4 (-21(1 - 2xf + 3)] + i (-x6 5 + (1 - x)6 6 )| 
where £ = 1 — 2x and the Gegenbauer Polynomials are given by, 



^ 2 (t) = *, ^ 2 (t) = i(3t 2 -i), ci'\t) = \(\e 



- 1 



C\ /2 (t) = - (3 - 30t 2 + 35t 4 ) 



1 



Cl'\t) = 3t , cY\t) = u -(5t* - 1) , CY\t) = £ (7t 3 - 3t) 



^3/2 



3/2/ 



(Bl) 



(B2) 



The coefficients are defined as 
61 



i + ^ + iSpX-yP-af, 



1 3 

63 = 5r] 3K - -r) 3K uJ3K + > & 4 = ViKhx 



5 5 (6) 



9(p*±p*)(l T 3af + 6af), p* = K + ^ , P- 

777- r 



(B3) 

"A" m A 

with m g being the mass of m n or since m n ~ is assumed. Since m q < m s , in our 
numerical estimations, we take p% = p K = p K ■ We display the values of decay constant, 
mass of strange quark, and relevant coefficients of distribution amplitudes for K meson at 
p = 1 GeV in Table El 
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TABLE V: The decay constant, mass of strange quark (in units of MeV) and coefficients of 
distribution amplitudes for K meson at n = 1 GeV. 



Sk 


m s 


af 


af 


P K 


r]3K 


^3K 


^3K 


160 


120 


0.06 


0.25 


0.076 


0.016 


-1.2 


1.6 



Similarly, the distribution amplitude for K* can be expressed as 33l. |34| 



<t>K*{x) 

4> P K »{x) 



f. 



K* 



2y/2N c 
fT 

Jk* 



6x{l-x) 1 + 3a^ + 3a 2 'C 2 3/2 (£) 



2^/2N~ c 



3e + 3a^C 2 1/2 (0 + 4C 2 3/2 (0 + 70C 3 T C 4 1/2 (O 



fl* 



+ (2 + ln(l -ar) + lnx) 



6^ 



+ + (^4 + |C 3 T ) (-203(1 - x) + 5£ 2 - l) 



-12a^ (1 - x) + 35+ (3f - 2 In (1 - x) - 2) 



(B4) 



The values of the decay constants and relevant coefficients of the distribution amplitudes 
for the K* meson are shown in Table I VI I 

TABLE VI: The decay constants (in units of MeV) and coefficients of distribution amplitudes for 
K* meson at \i = 1 GeV. 



fx* 


Jk* 


p> 


ajj 


«2~ 




5+ 


5- 


210 


170 


0.10 


0.09 


0.13 


0.024 


0.24 


-0.24 
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